Productions of hadrons, pentaquarks O"*^ and O*"*""*", and di-baryon (i7n)o+ in relativistic 
heavy ion collisions by a quark combination model 



Feng-Ian Shao,^'^ Qu-bing Xie,'^' ^ and Qun Wang"*'^ 

^Department of Physics, Shandong University, Jinan, Shandong 250100, People's Republic of China 
''Department of Physics, Qufu Normal University, 
Qufu, Shandong 273165, People's Republic of China 
^ China Center of Advanced Science and Technology (World Lab), 
P.O.Box 8730, Beijing 100080, People's Republic of China 
'^Institut fiir Theoretische Physik, J.W. Goethe-Universitat, D-60054 Frankfurt/Main, Germany 
^Department of Modern Physics, University of Science and Technology of China, 
Hefei, Anhui 230026, People's Republic of China 
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I _ - bination model. The mixed ratios for various hadrons and the transverse momentum spectra for 

' long-life hadrons are predicted and agree with recent RHIC data. The production rates for the 

, pentaquarks O"^, 0*^"^ and the di-baryon (r2f2)Q+ are estimated, neglecting the effect from the 

' transition amplitude for constituent quarks to form an exotic state. 
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I. INTRODUCTION 



ff^ ■ Recently a lot of data for hadron multiplicities have been published from the relativistic heavy ion collision (RHIC) 
J> ' at Brookhaven National Laboratory PjBIIiSIEIESB- One of the purpose of this experiment is to produce a 
OO , deconfined phase of quarks and gluons (QGP) under extreme conditions at high temperature and density by smashing 
two gold nuclei. Thus it is important to find reliable probes to judge if the QGP is really formed in the experiment, 
which is both theoretical and experimental challenges because one cannot directly detect the free quarks and gluons but 
their decay products. Hadron multiplicities and their correlations are observables encoding information on chemical 
\ properties of the medium generated from the heavy ion collisions. It is impressive that the statistical thermal models 
• provide a good description for the available multiplicity data at almost all energies in heavy ion collisions by only a 
few parameters llEIIllllllllllllllllllllllliallSIIi. While the ability to reproduce the mixed particle 
ratio from statistical models is not a proof that the emitting source is in thermal equilibrium 21], it is further 
evidence that such a state, which is one of the necessary conditions for creation of a QGP, has been created. The 
O , recombination or coalescence models [23, 123, lla, |2g, |23, |23, 123 are very successful in explaining the RHIC puzzles 
2 ' [S^l^ll^ll^l^ll^. e.g. an anomaly in the ratio of proton to pion which is unexpectedly high (reaching about 1) at 
' middle transverse momentum region. This implies the hadronization by quark recombination plays an important role 
^ ■ even in the hard regime characterized by high transverse momenta and provides a piece of evidence for decomfinement 
; (see, e.g. js^). mm 

. The earliest quark combination model (QCM) can be dated back to 1970s |33,|33|, which was proposed to describe 
d ' the multiparticle production or the hadronization in various reactions. Certainly the most popular hadronization 
models nowadays are the string model and the cluster model j39|,|40|,|^, but the great advantage of quark combination 
picture in describing the inclusive hadron production is its simplicity. The success of the quark combination model 
in almost all kinds of high energy collisions, e.g. electron-positron, hadron-hadron and nucleus-nucleus collisions, is 
partly due to the universal stochastic nature of fragmentation or hadronization. In this sense, the QCM resembles 
the thermal or the statistical model, but it encodes more microscopic information. 

We have developed a variant of the QCM based on a simple quark combination rule ^3iE3- Using our QCM, 
we have described most of multiplicity data for hadrons in electron-positron and proton- proton/anti- proton collisions 
p4, 45, 46, 47, 48, 49, ^J- Also we solved a difficulty facing other QCMs in describing the TASSO data for baryon- 
antibaryon correlation: they can be successfully explained by our QCM [isL IsH . Embedded in the event generator, 
our QCM can also reproduce most of the global properties of hadronic events such as momentum spectra in electron- 
positron collisions 52, 53, 54, 55, 56J. Encouraged by the success of the statistical and recombination models in heavy 
ion collisions, in this paper we try to extend our QCM to reproduce the recent RHIC data for hadron multiplicities. 
Especially we will predict the production rates of three exotic baryons: the pentaquarks 0"*" and 0*"'"+ and di-baryon 

The pentaquark G"*" is an exotic baryon made of five quarks uudds, which have been discussed in the context of 
quark models in the early days of QCD [STtlssf. In 1997, Diakonov et. al. [5^ predicted the mass and width of 0+ 
at about 1540 MeV and 15 MeV respectively, using the chiral soliton model. Recent works on the property of the 



2 



pentaquark in the chiral field model can be foun d, e.g. i n Ref. |6Ctl6lll62l |. Several experimental groups have reported 
the discovery of the pentaquark state O"*" ^3. 163. 165ll66ll67| . Another pentaquark 9*++ belongs to the 27-plet baryon 
with flavor content uuuds and possible spin-parity = (3/2)+ [gJ]. A search for the 0*++ pentaquark has been 
carried out by BaBar collaboration ^6§\ in the decay channel B+ ppK^. But the results in this experiment are 
preliminary yet to be confirmed. The last exotic state we are going to look at is the di-baryon (TMTW • For more than 
20 years the search for di-baryon has been another important attempt in hadronic phvsics |69l iTCt ItH l72l It^ l74j . One 
believes that if di-baryons do exist those with multi-strangeness must be ideal candidates to be observed because of 
their relatively long lifetime due to the stability with respect to strong decay. Recently the structure and properties of 
di-baryons with large strangeness are investigated in the chiral S'C/(3) model [zS Iz3j Izll which quite successfully 
reproduces several nuclear properties. They found that some six-quark states with high strangeness have considerable 
binding energy provided by the chiral quark coupling. Of particular interest is (rif2)g+ : it is a deeply bound state 
(binding energy is around 100 MeV; the mean-square root of the distance between two fi's is 0.84 fm). The mean 
lifetime of this di-baryon is as long as about twice that of f2, because it only decays weakly. All these interesting 
properties together with the electric charge Q — ~2 would make it easily identified in the experiments. Because of 
its large strangeness, (J7f2]n+ is not likely to be produced in proton-proton collisions. But one expects the enhanced 
strangeness production |l4l lisl in heavy ion collisions at RHIC energies, which can be a best place to study the 
production of the di-baryon (rjri)o+ |79ll80l|. 

The outline of this paper is as follows. In section ITU and ITTll we give a brief description of the quark combination 
model and basic relations among production weights of SUf{3) multiplets. In section Hvl we present particle ratios 
and the transverse momentum spectra of the pion, proton and kaon predicted from the quark combination model. 
The upper limits for production rates of 0+ and (1717)o+ are estimated in section IVl and IVTI We give a summary of 
results and conclude in section IVlTl 



II. THE QUARK COMBINATION MODEL 

All kinds of hadronization models demand themselves, consciously or not, satisfy rapidity or momentum correlation 
for quarks in the neighborhood of phase space. The essence of this correlation is its a gree ment with the fundamental 
requirement of QCD which uniquely determines the quark combination rule (QCR) M^. According to QCD, a qq 
may be in a color octet or a singlet which means a repulsive or an attractive interaction between them. The smaller 
the difference in rapidity for two quarks is, the longer the interaction time is. So there is enough time for a qq to be 
in a color singlet and form a meson. Similarly, a qq can be in a sextet or an anti-triplet. If its nearest neighbor is a g 
in rapidity, they form a baryon. If the neighbor is a q, because the attraction strength of the singlet is two times that 
of the anti-triplet, qq will win the competition to form a meson and leave a q alone. Our quark combination model 
is based on the above QCD requirements. When the transverse momentum of quarks are negligible, all q and q can 
always line up stochastically in rapidity. The QCR reads: 

1. Starting from the first parton {q or q) in the line; 

2. If the baryon number of the second in the line is of the different type from the first, i.e. the first two partons 
are either qq or qq, they combine into a meson and are removed from the line, go to point 1; Otherwise they are 
either qq or 7jq, go to the next point; 

3. Look at the third, if it is of the different type from the first, the first and third partons form a meson and are 
removed from the line, go to point 1; Otherwise the first three partons combine into a baryon or an anti-baryon 
and are removed from the line, go to point 1. 

Here is an example to show how the above QCR works; all partons line up in rapidity and combine to hadrons as 
follows 

91 929394 95969798 99910 9l 1 912 913 9l49l5 916 9l79l8 919920 

Miqiq^) B{q^q^q^) M{qeqj) B{qsqgqio) M{qj^j^qi2) Af ((713915) B((7i49i69i7) -B(9i89i992o) (1) 

We note that it is straightforward to define the QCR in 1-dimensional phase space, but it is more complicated to 
have it in two or three dimensional phase space, where one does not have an order or one has to define an order in a 
sophisticated way so that all quarks can combine to hadron in a particular sequence, see for example Ref. 

If the quarks and anti-quarks are stochastically arranged in rapidity, the probability distribution for N pairs of 
quarks and anti-quarks to combine into M mesons, B baryons and B anti-baryons is 

2iV(iV!)^(M + 2i?-l)! 
Xmb{N} (2iv)!Af!(B!)2 ^ ^nm+sb. (2) 
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TABLE L The weights of quark combination and those in terms of baryon multiplets. 



Flavor content 


Combination weight 


Production weight of baryon multiplets 


uuu 


1 




ddd 


1 


PlO 


sss 




PiqX^ 


uud 


3 


Pio + -Ps 


uus 


3A 


[pio + ps]A 


ddu 


3 


Pio + Ps 


dds 


3A 


bio + PsJA 


uss 


3A^ 


[pio + psJA^ 


dss 


3A2 


[pio + Ps]A^ 


uds 


6A 


[pio + 2(ps + Pi/)] A 



The average number of primarily produced mesons M{N) and baryons B{N) are given by 

{M{N)) = ^^A/XmbW, (3) 

M B 

{B{N)) = Y.Y.BXmb{N). (4) 

M B 

Approximately, for iV > 3, {M{N)) and {B{N)) can be well parameterized as linear functions of quark number N: 
{M{N)) ^aN + b and {B{N)) = (1 - a)7V/3 - b/3 where a = 0.66 and b = 0.56. But for iV < 3, obviously one has 
(M(iV)} = N and (B{N)) = 0. 

Having the number of mesons and baryons in an event, we can obtain the multiplicity of all primary hadrons from 
their production weights. The yield of the hadron hi can be written as 

(h,) = (M)Br(Mj h,) + Cb, {B)Bt{B, h,) + Y, {B)Bt(Bj ^ h,), (5) 

j j j 

where (M) and (B) are the average number of mesons and baryons respectively. Cmj j Cs^ a-nd C-g . are the normalized 
weights for the primary meson Mj , the primary baryon Bj , and the primary anti-baryon Bj , respectively. Obviously 
we have the property Cbj = Cg • ^'"(^i ~^ ^0 weighted branching ratio for hj to hi. 

The production weights Cm- and C^, for primary hadrons satisfy the SUf{S) symmetry with a strangeness sup- 
pression factor As for strange hadrons 83i, ^Mo £5] • To determine Cm, and Cb^, we need the ratio V/P of the vector 
{J^ = 1^) to the pseudoscalar (J^ = 0^) meson, and the ratio (3/2)+/(l/2)+ of the decouplet (J^ = (3/2)+) to the 
octet (J^ — (1/2)+) baryon. Assuming SU{6) symmetry, Anisovich et al. |33 gave V/P = 3. In this case the weights 
for all mesons except rj and rj' can be simply written as Ca/^ oc {2Ji + 1)X^* where Ji is the spin of M^, and is 
the number of strange quarks and/or anti-quarks in the meson. However there is a spin suppression effect associated 
with baryons. In the next section, we derive the two relations between the production weights for octet, decuplet and 
singlet baryons from the properties of hadronization in the quark combination scheme. 

III. 517/ (3) SYMMETRY AND FLAVOR CONSERVATION IN THE QUARK COMBINATION MODEL 

Hadronization is the soft process of the strong interaction and is independent of flavor, so the net flavor number 
remains constant during the process, which we call the property of flavor conservation. In the quark combination 
scheme, this means that the number of quarks of a certain flavor prior to hadronization equals to that of all primarily 
produced hadrons after it. The As-broken 5J7/(3) symmetry in hadron production means baryons or mesons in the 
same J^^ multiplet share an equal production rate up to a A^ factor. This SUf{3) symmetry has been supported by 
many experimen ts, p articularly by the fact that the observed As obtained from various mesons and baryons coincide 
with each other |83 |. This experimental fact is unexpected in the usual diquark model for baryon production and 
turns out to be in favor of the quark combination scheme |86| where a baryon is formed by the stochastic combination 
of three constituent quarks. 

From the above properties, we can obtain relations among the production weights for octet, decuplet and singlet 
baryons in the quark combination picture[46j|. As shown in Table I, all the flavor combinations are listed in the first 
column, with their combination weights in the second. In hadronization three quarks combine into a primary baryon 
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which satisfies the As-broken SUf{3) symmetry. As each flavor combination in the first column corresponds to a 
certain baryon belonging to several SUf(S) multiplets whose production weight is listed in the third column. The 
ground-state decuplet and octet baryons are denoted by 10 and 8, while the only excited baryon which we consider, 
the singlet A(1520), is denoted by 1'. Their corresponding weight are denoted by Pio, Ps, and respectively. Note 
that the two sets of weights must be associated with a common factor, we finally derive the following relations among 
the production weights for octet, decuplet and singlet baryons: 

Pio - Pv (6) 
Ps = 2Pio (7) 

They impose a global constraint on the production rates of all ground state baryons and the excited singlet baryon 
A(1520). We can understand the so-called spin suppression from the basic relations. For those decuplet and octet 
baryons which are primarily produced and have the same strangeness, the ratio of their production rate is 

(3/2)+ _ Pio _ 

Note that the ratio R measured in experiments is for the octet baryons which include decay products from the decuplet 
ones, we then obtain the following approximated value for Rexp when neglecting the production of excited baryons: 



(3/2)+ (ground) Pn 
(l/2) + (ground) ^ Pio + Ps 



Note that Rexp is much less than 2 from spin counting, which is called the spin suppression effect. 

We can see that there is an essential difference between the spin suppression factor R for the baryon and the 
multiplicity ratio of vector mesons to pseudo-scalar mesons V/P. The former is for the ratio of the (3/2)+ decuplet 
baryon to the (1/2)+ octet one. They belong to different SUf{3) multiplets. The basic relations impose a constraint 
upon their production weights. The latter is the ratio of vector to pseudo-scalar mesons, both of which belong to 
SUf(3) nonets. Flavor conservation and SUf{3) symmetry hold in each nonet and hence there is no restriction on 
their weights. 



IV. HADRON MULTIPLICITIES AND MOMENTUM SPECTRA 



In this section,we use the QCM to compute hadron multiplicities and their ratios in heavy ion collisions at RHIC 
energies. We will also calculate the transverse momentum spectra for pions, kaons and protons. Before we do that, 
we have to determine some input parameters of the QCM. The parameters which control the total multiplicity are the 
number of quarks and anti-quarks. In electron-positron and proton-antiproton collisions, the number of quarks is equal 
to that of anti-quarks, which means there are no excess baryons in contrast to anti-baryons. But for nucleus- nucleus 
collisions, we need two parameters, the number of quarks and that of anti-quarks, to account for the net baryon 
number even in central rapidity. These two parameters are determined by fitting the the total charged multiplicity 
data, {Nch)data = 4100 ± 210, in central Au-f Au collisions at 130 AGeV [l|, which corresponds to the total number 
of quarks and anti-quarks {Nq + Nq) = 7400 in the QCM. The quark number {Nq) and the anti-quark number {Nq) 
can be further determined by the ratio of antiproton to proton p/p = 0.7 3, 4]. Then we obtain (Nq) = 3920 and 
(Nq) = 3500, where we see that the net quark number is about 420. Another parameter is the strangeness suppression 
factor As which has also to be input from data. We find As = 0.5 is consistent with data ^/AT*" = 0.47 2]. At 200 
AGeV, we determine in the same way the total number of quarks and anti-quarks {Nq + Nq) = 8900 with the net 
quark number 360 and As — 0.6 by fitting the data [8ll |. 

Having determined the above parameters, we calculate the ratios of strange anti-baryons to their baryon counter- 
parts, A/A, S /S^,r2 /^^, in mid-rapidity of central Au+Au collisions and compare our results with STAR data 
I^Q. The above antibaryon to baryon ratios are mainly controlled by the net quark number from colliding nucleons. 
The results are listed in Table II. The data show that the ratios increase with strangeness of the baryons. The reason 
for this trend is that the quark pair production is more important than the baryon transport in mid-rapidity. We see 
that there is a good agreement between our model predictions and the datak, which means the quark combination 
mechanism can explain this behavior. In the QCM various ratios B/B are associated with a common multiplicative 
factor D which is given by the ratio K~^/K~, for example A/A = D{p/p), E/E = D{A/A) and 0/f2 = D{E/E). The 
calculated ratio K'^ /K~ and the IJ-factor in three compound ratios as shown above. The results are listed in Table 
III. We also calculate the ratios of singly-strange particles to non-strange particles and those of multi-strange particles 
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TABLE II: Anti-particle to particle ratios: our predictions and STAR data at 130 AGeV 0| and 200 AGeV fl. The second 
and fourth columns are STAR data at 130 and 200 AGeV respectively. The third and fifth columns are our results at 130 and 
200 AGeV respectively. 





STAR (130) 


QCM 


STAR (200) 


QCM 


p/p 


0.71±0.01±0.04 


0.71 


0.73 ± 0.05 


0.78 


A/A 


0.71±0.01±0.04 


0.79 


0.84 ± 0.05 


0.84 




0.83±0.04±0.05 


0.88 


0.94±0.08 


0.91 


n+/f2- 


0.95 ± 0.15 ± 0.05 


1.00 


1.03 ± 0.12 


1.00 



TABLE III: / K ratio compared with compound ratios of baryons at 130 AGeV 13). 





STAR 


QCM 


K+/K- 


1.092 ± 0.023 


1.122 


A/A 


0.98 ± 0.09 


1.10 


A/ A 


1.17 ±0.11 


1.11 




1.14 ±0.21 


1.15 



to singly-strange particles. The results and the experimental data |j, |5{ are given in Table IV and V. All above 
show that the mixed ratios for various hadrons agree with RHIC data. 

As a test of our QCM, we finally calculate the transverse momentum spectra for the pion, kaon and proton. Two 
steps are needed to get final hadron spectra. One is to obtain the spectra for initially produced hadrons which are those 
before their decay. We can do this through our Quark Combination Model. The next step is to let these initially 
produced hadron decay to final state hadrons which are observed in experiments. In our QCM, the longitudinal 
momentum is described by the constant distribution of rapidity, where we assume that quarks and anti-quarks are 
distributed in rapidity range y S [—2.8, 2.8] at 130 GeV and y £ [—3.0, 3.0] at 200 GeV with equal probability. This 
corresponds to the rapidity distribution of charged multiplicity at RHIC. Once the transverse momentum and rapidity 
for each quark are determined, its momentum is fixed. We can randomly line up all quarks in rapidity and let them 
combine into hadrons following our combination rule in section^ The momentum of a hadron is then the sum of that 
of its constituent quarks, while its energy is given by the mass-shell condition. The transverse momentum distribution 
of quark or anti-quark is extracted from the measured neutral pion spectra at 200 AGeV (87| , where we assume that 

they are identical because the observed ratios of anti-particles to particles, i.e. p/p, A/A, S /S^ and /K~ , are 
almost constant with pT This is what we get for quark or antiquark px distribution: f{pT) — {Pr^ + l)"'^ ". 

We assume that the combination occurs only among those constituent quarks which have the same azimuthal angle, 
similar to the strategy in Ref. |22l |. Note that all these hadrons after combination are primarily produced, in order 
to get the spectra comparable to data one has to let them decay in their center-of-mass system. The momenta of 
all decayed hadrons are then boosted back into laboratory frame and recorded. Here we make use of the approriate 
function of event generator JETSET. Finally we obtain the pr spectra for the pions, kaons and protons in final state 
which includes all contributions from decay. 

That we distinguish primarily produced hadrons from the final state ones marks the essential difference between 
our model and other recombination or coelescence models where one does not. Just 

to get the feeling of this substantial difference, one may take pions as an example: even in the range pt > 2GeV/c, 
about 70% of pions come from resonance decays. This is the reason why our quark px spectrum extracted from 
data decreases not exponentially at low transverse momenta but in a power law - it reflects the spectra of initially 
produced hadrons instead of final state ones. Another difference between our model and other combination models 
is that the quark pT distribution is given by one single function without distinguishing thermal quarks (low px) and 



TABLE IV: A group of ratios for strange hadrons compared with STAR data at 130 AGeV 





DATA 


QCM 




0.49 ± 0.05 ± 0.12 


0.47 


K+/Tr- 


0.161 ± 0.002 ± 0.024 


0.146 




0.146 ± 0.002 ± 0.022 


0.130 


(TT+ + n-)/h- 


(2.24 ± 0.69) X 10"'^ 


3.21x10"-' 
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TABLE V; A group of ratios for strange hadrons compared with STAR data at 200 AGeV (ii^. 





DATA 


QCM 




0.146 ± 0.024 


0.132 




0.09 ± 0.01 


0.089 


Q- /h- 


(1.20 ± 0.12) X 10"^ 


1.40 X 10"^ 


7r ^ / TT^ 


1.025 ± 0.006 ± 0.018 


1.008 


K-/K + 


0.95 ± 0.03 ± 0.03 


0.92 



TABLE VI: Upper limits for multiplicities of 9+, 9*++, and {nn)o+ at 130 AGeV and 200 AGeV. The second and third 
columns are predictions at 130 AGeV, while the fourth and fifth columns are predictions at 200 AGeV. Per event values are 
obtained by taking the full rapidity range as \y\ < 3.43. 





per rapidity (130) 


per event 


per rapidity (200) 


per event 


9+ 


1.16 


7.97 


1.30 


8.93 


9*++ 


1.86 


12.75 


2.08 


14.29 


fl- 


0.40 


2.74 


0.51 


3.50 


(m)o+ 


2.70x10"^ 


1.86x10"* 


3.67x10"^ 


2.52x10'* 




6.79 xlO"-"^ 


6.79x10"^ 


7.20x10"^ 


7.20x10"^ 



shower quarks (high pt)- In any combination models the role of gluons has been finally taken by quark- antiquark 
pairs, for example, shower quarks are developed after a QCD cascade. In our QCM, the number of constituent quarks 
already effectively includes those quark-antiquark pairs converting from gluons. All quarks in the whole pr range 
are assumed to combine following the same combination rule. This assumption is widely used in other combination 
models: combination among shower quarks is treated in the same footing as among shower quarks and thermal quarks. 
The difference just lies in the sources: shower quarks are products after QCD showering while thermal quarks are 
from thermal distribution. Except these differences, our model is basically the same as other combination models. 
This can be seen that the rapidity correlation principle in our model is also implied in other models: in Hwa- Yang's 
model, they assumed yi = y2 = Vh and {yi + y2)/2 — yh where yi and 2/2 are quark rapidities, and yu is the hadron 
one [23; While in GKL model |2J, they used even more simple assumption yi = 2/2 = J/ft, = 0. 

Instead of giving the prediction for pT spectra, we give the spectra of the scaling variable z — pr/K where the 
scaling factor K is defined in Ref.|22]- The result for pions, kaons and protons at 130 and 200 AGeV are shown in 
Fig. H1I3|I . We see that the agreement between our predictions and data is quite satisfactory and we have verified 
the perfect scaling behavior in our model. The observed enhancement of protons and antiprotons at intermediate 
transverse momenta can be also naturally explained because we have got right predictions for protons. 

V. MULTIPLICITIES OF 9+ AND 9*++ 

In this section we will give an estimate for multiplicities of 0+ and 0*++ at 130 AGeV and 200 AGeV. Several 
groups have already predicted the yields of 0+ in central Au+Au collisions by the statistical and coalescence models 
[89, 90, 91]. 

In the QCM it is quite easy to estimate the production rate of 9+ and 9*++. The probabilities P{uudds) and 
P{uuuds) for five quarks uudds and uuuds respectively to come together in rapidity can be easily obtained by the 
QCM. The overlapping of phase space (here rapidity) has been encoded in the probability P. Then the probabilities 
to form the pentaquark 9'^ and 9*"'""'" can be estimated by the spin counting: a system composed of five quarks, each 
of which has spin-1/2, has totally 2^ spin states composed of one spin-5/2, four spin-3/2 and five spin-1/2 states. 
Assuming that all spin-1/2 states form 9+, the production rate for 9+ is then P(9+) = ^P{uudds). If all spin-3/2 
states go to 9*"'""'", the production rate of 9*+"'" is then P(9*"'""'") = ^P{uuuds). Note that the current spin counting 
method only provides a kind of upper limits for the yields. In the real calculation one must evaluate the transition 
probability which involves overlapping of the wave functions of pentaquark states and their quark constituents. Here 
we simplify the problem by assuming the transition amplitude as unity. Also we neglect the contribution from the 
hadron-hadron rescattering in hadronic phase and assume the dominant source is from quark combination. Our 
estimates for 9+ and 9*+^^ are shown in Table VI. In summary of our results, at 130 AGeV there are about 1.16 
9+ and 1.86 9*''"+ per rapidity produced in mid-rapidity in central Au-|-Au collisions, while at 200 AGeV, the yields 
of 9+ and 9*++ are 1.3 and 2.08 per rapidity respectively. Our predictions agree with those given by the statistical 



7 



models [s^ [oOl- For comparison, the yield of O"*" in proton- proton collision is predicted to be about 10 ^ — 10 

mill. 

VI. MULTIPLICITY OF {m)o+ 

We have calculated the ratios of strange anti-baryons to their baryon counterparts and that of to negative charged 
particles. The results are found to be consistent with available data, see Table IV-V. Based on the result for fi, we 
will predict in this section multiplicities of the di-baryon (riri)Q+ in Au+Au collisions at 130 and 200 AGeV. 

First we estimate the yield of ft {fl^ + fl^). Assuming the ratio of fl to h~ is constant for all rapidities, we get 
the yield of ft at about 5 and 6 per event from the experimental data h~ ~ \/2N'^^*'°' w 2050 at 130 AGeV and 
/i^ w 1/2 A*"^^*" w 2450 at 200 AGeV lU respectively. In the quark combination picture, the production weight for 
is proportional to A^. We find that the multiplicities increase strongly with growing A^, from 0.76 at As = 0.3 to 
6.00 at As = 0.7. Our predictions for the Vl^ yields are listed in Table VI. One sees that the yield of Vl~ per event 
is about 2.74. Given the result for il^, we are now in a position to estimate the multiplicity of (r2r2)o+. First we 
get the probability P{ssssss) for six s-quarks to come together under the rule of our QCM and then we apply the 
spin-counting to estimate the yield of (r2r2)o+ . There are totally 2^ spin states for a six-quarks system. In terms of 
total spin, these are one spin-3, five spin-2, nine spin-1 and five spin-0 states. We assume that all spin-0 states go to 
(riri)o+ . The yield of (riri)o+ is then ^ times that of six s quarks. The results are shown in Table VI. Same as the 
results for pentaquarks, the yields of (f2f7)o+ in Table VI are only regarded as a kind of upper limits because we have 
neglected the effect of the transition amplitude for 6 s-quarks to form a di-omega. One sees that the multiplicities of 
(riri)o+ are about 1.86 x 10~^ at 130 AGeV and 2.52 x 10~^ at 200 AGeV in an event respectively. Our predictions 
agree in magnitude with that of the coalescence models [t^ and are within the scope of present RHIC experiments. 
Normally one can detect (ilf2)o+ through its decay product. Because the binding energy of (1717)o+ is relatively large 
which prevent it from strong decay, (f2f2)o+ mainly weakly decays into one plus the decay products of (three or 
more body decay) or into Vl and S (two body decay). If it is produced in the experiments, one could detect it through 
the above decay channels. For other experimental traces of dibaryons, see, for example, Ref. [o^. l95|. 

VII. SUMMARY 

In this paper, we extend our quark combination model, which is very successful in describing the hadron production 
in electron-positron and proton-proton(anti-proton) collisions, to reproduce the multiplicity data in heavy ion collisions 
at RHIC energies 130 and 200 AGeV. The model can describe available data for transverse momentum spectra for 
pions, kaons and protons, especially the px scaling behavior can be well reproduced within our simple model. It can 
explain the anormaly of the proton to pion ratio at intermediate transverse momenta. A variety of ratios for anti- 
baryons to baryons, single-strange hadrons to non-strange ones, and double-strange hadrons to single-strange ones 
can be all reproduced in our QCM. With the approximation that the transition probability for constituent quarks to 
form an exotic state is taken as unity if they come togather in rapidity according to our combination rule, we finally 
estimate the yields of 9+, 0*++ and (f2fi)o+ in Au-f Au collisions at 130 and 200 AGeV. These yields can be regarded 
as a kind of upper limits. The multiplicities of O"*" and 0*^*"+ are estimated to be of order 1 and 2 per rapidity in 
mid-rapidity in a central event respectively, while the rate of (Sm)o+ is of the magnitude 10~^ per rapidity. The 
multiplicites of O"^ and O*"*""*" are found to be almost independent of the strangeness suppression factor As, while that 
of (f7r2)o+ increases strongly with growing As. 
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